Search for: All records

A new version of the US National Science Foundation National Center forAtmospheric Research (NSF NCAR) thermosphere-ionosphere-electrodynamicsgeneral circulation model (TIEGCM) has been developed and released. Thispaper describes the changes and improvements of the new version 3.0since its last major release (2.0) in 2016. These include: 1) increasingthe model resolution in both the horizontal and vertical dimensions, aswell as the ionospheric dynamo solver; 2) upward extension of the modelupper boundary to enable more accurate simulations of the topsideionosphere and neutral density in the lower exosphere; 3) improvedparameterization for thermal electron heating rate; 4) resolvingtransport of minor species N(2D); 5) treating helium as a major species;6) parameterization for additional physical processes, such as SAPS andelectrojet turbulent heating; 7) including parallel ion drag in theneutral momentum equation; 8) nudging of prognostic fields near thelower boundary from external data; 9) modification to the NO reactionrate and auroral heating rate; 10) outputs of diagnostic analysis termsof the equations; 11) new functionalities enabling model simulations ofcertain recurrent phenomena, such as solar flares and eclipse. Wepresent examples of the model validation during a moderate storm andcompare simulation results by turning on/off new functionalities todemonstrate the related new model capabilities. Furthermore, the modelis upgraded to comply with the new computer software environment at NSFNCAR for easy installation and run setup and with new visualizationtools. Finally, the model limitations and future development plans arediscussed. 

An expression of Rayleigh‐Taylor (R‐T) instability growth rate based on the field‐line integrated theory is newly established. This expression can be directly utilized in ionosphere models with magnetic flux tube structure based on Modified Apex Coordinates. The R‐T instability growth rates are calculated using the thermospheric and ionospheric conditions based on the coupled Whole Atmosphere Model and Ionosphere Plasmasphere Electrodynamic model (WAM‐IPE). The parameters used in this calculation include the field‐line integrated conductivities and currents, which consider the Quasi‐Dipole Coordinates and the modifications to the equations of electrodynamics. Detailed description of the new formulas and comprehensive analyses of diurnal, longitudinal, and seasonal variations of the R‐T instability growth rate are carried out. The dependencies of growth rates on pre‐reversal enhancement (PRE) vertical drifts and solar activity are also examined. The results show that pronounced R‐T growth rates are captured between 18 and 22 local time (LT) when strong PRE occurs in the equatorial ionosphere. The simulated R‐T growth rate increases with increasing solar activity levels and demonstrates strong correlations with the angle between the sunset terminator and the geomagnetic field line. These results are consistent with plasma irregularity occurrence rates shown in various satellite observations, suggesting that the newly developed R‐T growth rate calculation can effectively capture the probability of irregularities by considering the changes along magnetic flux‐tubes in the ionosphere. Since the WAM‐IPE is running in operation at National Oceanic and Atmospheric Administration Space Weather Prediction Center, the new calculations can be potentially implemented in the near future to provide forecasted information of the R‐T growth rate. 

Abstract We report the first lidar observations of regular occurrence of mid‐latitude thermosphere‐ionosphere Na (TINa) layers over Boulder (40.13°N, 105.24°W), Colorado. Detection of tenuous Na layers (∼0.1–1 cm⁻³from 150 to 130 km) was enabled by high‐sensitivity Na Doppler lidar. TINa layers occur regularly in various months and years, descending from ∼125 km after dusk and from ∼150 km before dawn. The downward‐progression phase speeds are ∼3 m/s above 120 km and ∼1 m/s below 115 km, consistent with semidiurnal tidal phase speeds. One or more layers sometimes occur across local midnight. Elevated volume mixing ratios above the turning point (∼105–110 km) of Na density slope suggest in situ production of the dawn/dusk layers via neutralization of converged Na⁺layers. Vertical drift velocity of TINa⁺calculated with the Ionospheric Connection Explorer Hough Mode Extension tidal winds shows convergent ion flow phases aligned well with TINa, supporting this formation hypothesis.